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Abstract 



^ , Starting from an effective Skyrme interaction we present a method to take 

into account the couphng between one- and two-phonon terms in the wave 
functions of excited states. The approach is a development of a finite rank 
separable approximation for the quasiparticle RPA calculations proposed in 
our previous work. The influence of the phonon-phonon coupling on energies 
and transition probabilities for the low-lying quadrupole and octupole states 
in the neutron-rich Sn isotopes is studied. 
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I. INTRODUCTION 



The experimental and theoretical studies of properties of the excited states in nuclei far 
from the /9-stability line are presently the object of very intensive activity. The random 
phase approximation (RPA) [1-4] is a well-known and successful way to treat nuclear vi- 
brational excitations. Using the Gogny's [5] or Skyrme-type [6] effective nucleon-nucleon 
interactions the most consistent models can describe the ground states in the framework of 
the Hartree-Fock (HF) and Hartree-Fock-Bogoliubov (HFB) approximations and the excited 
states within the RPA and quasiparticle RPA (QRPA). Such models are quite successful not 
only to reproduce the nuclear ground state properties [7,8], but also to describe the main 
features of nuclear excitations in closed-shell [9,10] and open-shell nuclei [11-14]. In the 
latter case the pairing correlations are very important. 

Due to the anharmonicity of vibrations there is a coupling between one-phonon and more 
complex states [2,4] and the complexity of calculations beyond standard RPA or QRPA in- 
creases rapidly with the size of the configuration space, so one has to work within limited 
spaces. Making use of separable forces one can perform calculations of nuclear charac- 
teristics in very large configuration spaces since there is no need to diagonalize matrices 
whose dimensions grow with the size of configuration space. For example, the well-known 
quasiparticle-phonon model (QPM) [4] can do very detailed predictions for nuclei away from 
closed shells [15], but it is very difficult to extrapolate the phenomenalogical parameters of 
the nuclear hamiltonian to new regions of nuclei. 

That is why a finite rank approximation for the particle-hole (p-h) interaction result- 
ing from the Skyrme forces has been suggested in our previous work [16]. Thus, the self- 
consistent mean field can be calculated with the original Skyrme interaction whereas the 
RPA solutions would be obtained with the finite rank approximation to the p-h matrix 
elements. It was found that the finite rank approximation can reproduce reasonably well 
the dipole and quadrupole strength distributions in Ar isotopes. Alternative schemes to 
factorize the p-h interaction were considered in [17-19]. 
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Recently, the finite rank approximation for p-h interactions of Skyrme type has been 
generahzed to take into account the pairing correlations [20]. The QRPA was used to 
describe characteristics of the low-lying 2+ and 3~ states and giant resonances in nuclei with 
very different mass numbers [20,21]. It was found that there is room for the phonon-phonon 
couphng effects in many cases. The first calculation to estimate this effect has been done 
for ii^Sn in [22]. 

In the present work, we extend our approach to take into account the coupling between 
the one- and two-phonon terms in the wave functions of excited states. As an application 
of the method we present results for low-lying 2+ and 3~ states in neutron-rich Sn isotopes 
and compare them with recent experimental data [23] and other calculations [24-26]. 

This paper is organized as follows: in Sec. II we sketch our method allowing to consider 
effects of the phonon-phonon coupling. In Sec. Ill we discuss details of calculations and 
show how this approach can be applied to treat the low-lying collective states. Results of 
calculations for properties of the quadrupole and octupole states in ^^^~^^^Sn isotopes are 
given in Sec. IV. Conclusions are drawn in Section V. 

II. METHOD OF CALCULATIONS 

A. The model hamiltonian and QRPA 

We start from the effective Skyrme interaction [6] and use the notation of Ref. [27] 
containing explicit density dependence and all spin-exchange terms. The single-particle 
spectrum is calculated within the HF method. The continuous part of the single-particle 
spectrum is discretized by diagonalizing the HF hamiltonian on a harmonic oscillator basis 
[28] . The p-h residual interaction Vres corresponding to the Skyrme force and including both 
direct and exchange terms can be obtained as the second derivative of the energy density 
functional with respect to the density [29]. Following our previous papers [16] we simplify 
Vres by approximating it by its Landau-Migdal form. For Skyrme interactions all Landau 
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parameters Fi,Gi, F^, G\ with I > 1 are zero. Here, we keep only the I — terms in Vres and 
in the coordinate representation one can write it in the following form: 

S{r,-r,) (1) 

where Uj and Tj are the spin and isospin operators, and A^o = ^kprn* /n'^h'^ with kp 
and m* standing for the Fermi momentum and nucleon effective mass. The expressions for 
Fq, Go, -Fq, G'q in terms of the Skyrme force parameters can be found in Ref. [27]. Because 
of the density dependence of the interaction the Landau parameters of Eq.(l) are functions 
of the coordinate r. 

In what follows we use the second quantized representation and Vres can be written as: 

Vres = ^ ^1234 : 0^^0^0403 : (2) 
^ 1234 

where af (oi) is the particle creation (annihilation) operator and 1 denotes the quantum 
numbers (nilijimi), 

Vi234 = y"0t(ri)02(r2)K-e5(ri,r2)03(ri)04(r2)dridr2. (3) 

After integrating over the angular variables one needs to calculate the radial integrals. 
As it is shown in [16,20] the radial integrals can be calculated accurately by choosing a large 
enough cutoff radius R and using a A^-point integration Gauss formula with abscissas 
and weights Wf.. Thus, the two-body matrix element is a sum of A^" separable terms, i.e., the 
residual interaction takes the form of a rank separable interaction. 

We employ a hamiltonian including an average HF field, pairing interactions, the isoscalar 
and isovector particle-hole (p-h) residual forces in a finite rank separable form [20] : 

^ = E ( (E, - Xr)a]^a,^ - \vl'^ : (r)Po (r) : ) + Ves, (4) 

r \jm ^ J 

where 

n+M= E' (-irXm<— (5) 

jm 
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No' 



Fo(ri) + Go(ri)ai ■ + (Po(ri) + Go(ri)cri ■ cr2)ri ■ T2 



We sum over the proton (p) and neutron(n) indexes and the notation {r = {njp)} is 
used. A change r — r means a change p ■(-^ n. The single-particle states are specified by 
the quantum numbers (jm), Ej are the single-particle energies, the chemical potentials. 
Vj>^^ is the interaction strength in the particle-particle channel. The hamiltonian (4) has the 
same form as the QPM hamiltonian with N separable terms [4,30], but the single-particle 
spectrum and parameters of the p-h residual interaction are calculated making use of the 
Skyrme forces. 

In what follows we work in the quasiparticle representation defined by the canonical 
Bogoliubov transformation: 

Qjm = uja^^ + {-ly-'^Vjaj-rn- (6) 

The hamiltonian (4) can be represented in terms of bifermion quasiparticle operators 
and their conjugates [4]: 

mm 

^^(ii'; V) = ^{jmj'm I A//)at^Q;J (8) 

mm 

We introduce the phonon creation operators 

QU = ^ E A^jf; Xn) - {-ir-^Y>^ Aijj'; A - /.)) . (9) 
jj' 

where the index A denotes total angular momentum and /i is its z-projection in the 
laboratory system. One assumes that the ground state is the QRPA phonon vacuum | 0) , 
i-e. Qxni I 0) =0. Wc define the excited states for this approximation by Qx^i \ 0). The 
quasiparticle energies (sj), the chemical potentials (A^), the energy gap and the coefficients 
u,v of the Bogoliubov transformations (6) are determined from the ECS equations with the 
single-particle spectrum that is calculated within the HF method with the effective Skyrme 
interaction. Making use of the linearized equation-of-motion approach [1] : 

0), (10) 
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SQxni, Q 



Xni 



with the normahzation condition: 



(0| [Q>.,i,QX,A |0) =5,/, 



(11) 



one can get the QRPA equations [3,4] : 



A E 
B -A 















= W 




I 









(12) 



In QRPA problems there appear two types of interaction matrix elements, the 
matrix related to forward-going graphs and the B^'f'., -, i, ,m matrix related to backward- 
going graphs. Solutions of this set of linear equations yield the eigen-energies and the 
amphtudes X, Y of the excited states. The dimension of the matrices A, B is the space 
size of the two-quasiparticle configurations. For our case expressions for A, B and X, Y are 
given in [20]. Using the finite rank approximation we need to invert a matrix of dimension 
X independently of the configuration space size [16,20]. Therefore, this approach 
enables one to reduce remarkably the dimensions of the matrices that must be inverted to 
perform structure calculations in very large configuration spaces. 

B. Phonon-phonon coupling 

Our calculations [20] show that, for the normal parity states one can neglect the spin- 
multipole terms of the p-h residual interaction (1). Using the completeness and orthogonality 
conditions for the phonon operators one can express bifermion operators A'^ [j j ; Xji) and 
A(jj'; A/i) through the phonon ones and the initial hamiltonian (4) can be rewritten in terms 
of quasiparticle and phonon operators in the following form: 



H = ho + hnn + h. 



QB 



3 '^jm '^J> 



(13) 
(14) 



jm 
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(15) 



hQB = E E^rjj, (r) {{-)'-'QU + Qx-,i) B{jj';X - /.) + h.c. (16) 

The coefficients W, T oi the hamiltonian (13) are sums of N combinations of phonon 
amphtudes, the Landau parameters, the reduced matrix element of the spherical harmonics 
and radial parts of the HF single-particle wave function (see Appendix A). It is worth to 
point out that the term Hqe is responsible for the mixing of the configurations and, therefore, 
for the description of many characteristics of the excited states of even-even nuclei [4] . 

To take into account the mixing of the configurations in the simplest case one can write 
the wave functions of excited states as: 



^.(A/^) = {E^i('^^)Q 

i 

with the normalization condition: 



AiiiA2i2 



(17) 



i AiiiA2i2 



Using the variational principle in the form: 



(18) 



S ((*4A/.)|i/ I ^.(A/.)) - E,((*.(A/x)|^.(A/.)) - 1)) = 0, 



(19) 



one obtains a set of linear equations for the unknown amplitudes Ri{Jv) and P\^il{Ji^)'- 



AljlA2«2 



(20) 



E-^ 



(21) 



U'^lll^Ji) is the matrix element coupling one- and two-phonon configurations [4,31]: 



(22) 
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The expression of U^lHiJi) is given in Appendix B. The number of hnear equations (20), 
(21) equals the number of one- and two-phonon configurations included in the wave function 
(17). 

The energies of excited states E^, are solutions of the secular equation 



F{E^) = det 



1 



Alii 



{uxi — E„)5iii — - ^ 



E, 



0, 



(23) 



Ain,A2i2 +^A2i2 

where the rank of the determinant equals the number of the one-phonon configurations. 
Using Eqs.(20), (21) and the normalization condition (18), one can find the amplitudes 
R,{Jv) ^^dP^^liJv). 

It is necessary to point out that the equations derived above have the same form as the 
basic QPM equations [4,31], but the single-particle spectrum and the p-h residual interaction 
are determined making use of the Skyrme interactions. 



III. DETAILS OF CALCULATIONS 

We apply the present approach to study characteristics of the low-lying vibrational states 
in the neutron-rich Sn isotopes. In this paper we use the parametrization SLy4 [32] of the 
Skyrme interaction. This parametrization was proposed to describe isotopic properties of 
nuclei from the /3-stability line to the drip lines. Spherical symmetry is assumed for the HF 
ground states. 

The pairing constants are fixed to reproduce the odd-even mass difference of neigh- 
boring nuclei. It is well known [11,12] that the constant gap approximation leads to an 
overestimating of occupation probabilities for subshells that are far from the Fermi level and 
it is necessary to introduce a cut-off in the single-particle space. Above this cut-off subshells 
don't participate in the pairing effect. In our calculations we choose the BCS subspace to 
include all subshells lying below 5 MeV. 

In order to perform QRPA calculations, the single-particle continuum is discretized [28] 
by diagonalizing the HF hamiltonian on a basis of twelve harmonic oscillator shells and 
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cutting off tlie single-particle spectra at the energy of 100 MeV. This is sufficient to exhaust 
practically all the energy-weighted sum rule. 

The Landau parameters Fq^Gq, Fq, G'q expressed in terms of the Skyrme force parameters 
[27] depend on kp- As it is pointed out in our previous works [16,20] one needs to adopt 
some effective value for kp to give an accurate representation of the original p-h Skyrme 
interaction. For the present calculations we use the nuclear matter value for kp. 

Our previous investigations [20] enable us to conclude that N—Ab for the rank of our 
separable approximation is enough for multipolarities A < 3 in nuclei with ^4 < 208. In- 
creasing N, for example, up to A^=60 in ^osp^^ changes results for energies and transition 
probabilities not more than by 1%. Our calculations show that, for the natural parity states 
one can neglect the spin-multipole interactions and this reduces by a factor 2 the total 
matrix dimension, i.e., the matrix dimensions never exceed 2N x 2N independently of the 
configuration space size [16,20]. 

The two-phonon configurations of the wave function (17) are constructed from natural 
parity phonons with multipolarities A = 2,3,4,5. All one-phonon configurations with en- 
ergies below 8 MeV for i24-i30,i34gj^ MeV for ^^^Sn are included in the the wave 
function (17). The cut-off in the space of the two-phonon configurations is 21 MeV. An 
extension of the space for one- and two-phonon configurations does not change results for 
energies and transition probabilities practically. 

IV. RESULTS OF CALCULATIONS 

As an application of the method we investigate effects of the phonon-phonon coupling 
on energies and transition probabihties to 2f and 3^ states in ^^^~^^^Sn. 

Results of our calculations for the 2f energies and transition probabihties B(E2) are 
compared with experimental data [23,33] in Table I. Columns "QRPA" and "2PH" give 
values calculated within the QRPA and taking into account the phonon-phonon coupling, 
respectively. 
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As it is seen from Table I there is a remarkable increase of the 2^" energy and B{E2 f) 
in ^^^Sn in comparison with those in -"^^^'^^^Sn. Such a behaviour of B{E2 |) is related 
with the proportion between the QRPA amplitudes for neutrons and protons in Sn isotopes. 
The neutron amplitudes are dominant in all Sn isotopes and the contribution of the main 
neutron configuration {lhn/2j l/'ii/2} increases from 81.2% in -"^^^Sn to 92.8% in -"^^^Sn when 
neutrons fill the subshell l/tii/2- At the same time the contribution of the main proton 
configuration {2^5/2, lg'9/2} is decreasing from 9.3% in ^^^Sn to 3.9% in ^^°Sn. The closure of 
the neutron subshell l/iii/2 in ^^^Sn leads to the vanishing of the neutron paring. The energy 
of the first neutron two-quasiparticle pole {2/7/2, l/iii/2} in ^^^Sn is greater than energies 
of the first poles in i30,i34gj^ ^j^g contribution of the {2/7/2, l/iii/2} configuration in the 
doubly magic ^^^Sn is about 61%. Furthermore, the first pole in ^^^Sn is closer to the proton 
poles. This means that the contribution of the proton two-quasiparticle configurations is 
greater than those in the neighbouring isotopes and as a result the main proton configuration 
{2(^5/2, ^99/2} in ^^^Sn exhausts about 33%. In ^^^Sn the leading contribution (about 99%) 
comes from the neutron configuration {2/7/2,2/7/2} and as a result the B{E2) value is 
reduced. Such a behaviour of the 2^ energies and B{E2) values in the neutron-rich Sn 
isotopes refiects the shell structure in this region . It is worth to mention that the first 
prediction of the anomalous behaviour of 2+ excitations around ^"^^Sn based on the QRPA 
calculations with a separable quadrupole- plus-pairing hamiltonian has been done in [24]. 
Other QRPA calculations with Skyrme [13] and Gogny [26] forces give similar results for Sn 
isotopes. 

One can see from Table I that the inclusion of the two-phonon terms results in a decrease 
of the energies and a reduction of transition probabilities. Note that the effect of the two- 
phonon configurations is important for the energies and this effect becomes weak in ^^^Sn. 
There is some overestimate of the energies for the QRPA calculations and taking into account 
of the two-phonon terms improves the description of the 2^ energies . The reduction of the 
B[E2) vahics is small in most cases due to the crucial contribution of the one-phonon 
configuration in the wave function structure. 
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Results of our calculations for the 3i energies and the transition probabilities B{E3) 
compared to experimental data [34] are shown in Table II. As for the quadrupole excitations 
the influence of coupling between one- and two-phonon terms in the wave functions of the 
3i states leads to the decrease of the energies and the reduction of transition probabilities. 
In spite of the fact that the 3i states have strong collectivity and many two-quasiparticle 
configurations give a contribution in the QRPA wave functions in Sn isotopes the phonon- 
phonon coupling is not very strong in this case. Our calculation shows that the main reason 
is the smallness of the matrix elements coupling the one-phonon configuration {3i} and the 

_i_ _ 2+ _ _ 

two-phonon configuration {25^;3i } {U J. (S^ )). As a result the decrease of the 3^ energies 
is about 10%. In the present paper we neglect the p-p channel that can be important for 
collective phonons and can reduce the collectivity of states [4,35]. This can give an additional 
lowering of energies and transition probabilities , but this is not the case for ^^^Sn. It is 
worth to mention that experimental data for 3i states in the neutron-rich Sn isotopes are 
very scarce. 

An additional information about the structure of the first 2+, 3~ states can be extracted 
by looking at the ratio of the multipole transition matrix elements Mn/Mp that depends 
on the relative contributions of the proton and neutron configurations. In the framework of 
the collective model for isoscalar excitations this ratio is equal to Mn/Mp = N/Z and any 
deviation from this value can indicate an isovector character of the state. The M^/Mp ratio 
can be determined experimentally by using different external probes [36-38] . Our calculated 
values for the Mn/Mp ratios for the 2i and 3j~ states arc shown in Table III. The calculated 
Mn/Mp ratios are rather close to N/Z except 2^ in ^^^Sn. It is worth to note that the 
deviation of the ratio for 2^ state in ^^^Sn correlates with the increase of the contribution 
of the proton two-quasiparticle configurations. 
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V. CONCLUSIONS 



A finite rank separable approximation for the QRPA calculations with Skyrme interac- 
tions that was proposed in our previous work is extended to take into account the coupling 
between one- and two-phonon terms in the wave functions of excited states. The suggested 
approach enables one to reduce remarkably the dimensions of the matrices that must be 
diagonalized to perform structure calculations in very large configuration spaces. As an 
application of the method we have studied the behavior of the energies and transition prob- 
abihties to 2^ and 3^ states in ^^^~^^^Sn. The inclusion of the two-phonon configurations 
results in a decrease of the energies and a reduction of transition probabihties. It is shown 
that there is some overestimate of 2i energies for the QRPA calculations and the effect of 
the two-phonon configurations is important, but this effect decreases in ^^^Sn. The inclu- 
sion of the two-phonon terms does not change the effect of a remarkable increase of the 
QRPA value of B{E2;0~^ 2^) for the doubly-closed shell nucleus ^^^Sn in comparison 
with its neighbors. A systematical study of the infiuence of the two-phonon terms taking 
into account the p-p channel on properties of the low-lying states is now in progress. 
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APPENDIX A 



The coefficients of the hamiltonian (13) are given by the following expressions: 




(Al) 



12 



(A2) 



where 



2{2\ + l] 



3J 



(+) 



2 ' 



In the above expressions, f .., denotes the single-particle radial matrix elements [20]: 



fhS = %i(^fc)«i2(^fc)^^(jl||>A||j2), 

where Uj^{rk) is the radial part of the HF single-particle wavefunction at the abscissas of 
the A^-point integration Gauss formula r^. k'q^'''^ and Ki^'''^ are defined by the Landau 
parameters as 



/ {M,k) \ 



JM,k) 



= -N, 



_i Rwk 



2rl 



Foir-k) 



APPENDIX B 



The matrix elements U^^-^(Ji) have the following form: 




V j2i3 isil jail i2i3 ^ 



/-y^A2j2"y^Aj I "V^Ai ■V^A2i2\ I 
\ jsjl 3233 3233 3331 ) 

\ 
/ 



' ' 3233 3331 ~^ 3331 3233 



(Bl) 
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TABLES 



TABLE I. Energies and B(E2)-values for up-transitions to the first 2+ states 



Nucleus Energy B(E2t) 

(MeV) (e^b^) 
Exp. Theory Exp. Theory 







QRPA 


2PH 




QRPA 


2PH 




1.13 


1.92 


1.03 


0.1660±0.0040 


0.177 


0.151 


i26Sn 


1.14 


1.96 


1.30 


0.10±0.03 


0.149 


0.133 


i28Sn 


1.17 


2.08 


1.48 


0.073±0.006 


0.111 


0.100 


Sn 


1.22 


2.37 


1.73 


0.023±0.005 


0.064 


0.058 


i32Sn 


4.04 


4.47 


4.03 


0.14±0.06 


0.136 


0.129 




0.73 


1.65 


1.34 


0.029±0.006 


0.016 


0.015 
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TABLE II. Energies and B(E3)-values for up-transitions to the first 3 states 



Nucleus 




Energy 




B(E3T) 








(MeV) 




(e^b^) 






Exp. 


Theory Exp. 


Theory 








QRPA 


2PH 


QRPA 


2PH 




2.60 


3.64 


3.25 0.073±0.010 


0.208 


0.196 


i26Sn 


2.72 


4.16 


3.76 


0.191 


0.176 






4.66 


4.22 


0.181 


0.161 






5.17 


4.75 


0.183 


0.159 


i32Sn 


4.35 


5.66 


5.36 


0.202 


0.191 


on 




5.01 


4.51 


U.IZO 


nil! 




TABLE III. (Mn/Mp) 


/(N/Z) ratios for the first 2+, 3^ 


" states 




State 


i24Sn 


i26Sn 


i28Sn i30sn 


i32Sn 




2^ 


0.99 


0.99 


0.98 0.97 


0.81 


1.44 


3r 


0.94 


0.92 


0.89 0.86 


0.83 


0.91 
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